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^ ; Abstract 

(N 

Hwang in [Phys. Rev. E 80, 042103 (2009)] suggested that the two transitions of the six-state 

Q ■ clock model on the square lattice are not of the Kosterlitz-Thouless type. Here we show from 

^ . 

^ . simulations that at the upper transition, the helicity modulus does make a discontinuous jump to 



zero. This gives strong evidence for a Kosterlitz-Thouless transition. 
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Recently, Hwang [l| has examined the six-state clock model by the Fisher-zero approach 
and suggested that its phase transitions are not of the Kosterlitz-Thouless (KT) type, in 
contrary to earlier theoretical analyses |2|, y]. This raises questions of earlier numerical works 
where scaling indices have been measured and found to be in agreement with the earlier the- 
oretical KT predictions (see, for example, Ref. |J]). In spite of the theoretical and numerical 
supports for the KT scenario in the six-state clock model, one might perhaps still argue that 
the actual transitions can nevertheless be of a standard continuous type due to the follow- 
ing reasons. First, the theoretical predictions usually relate to the Villain approximation 
which is not an exact representation of the clock model. Second, the agreements with the 

KT scaling indices only provide a necessary condition but maybe not a sufficient one for 

I, 
ruling out a standard continuous transition. Third, it has been claimed in Ref. [3] that the 

six-state clock model does not exhibit a discontinuous jump in the helicity modulus (T), a 
key feature of the KT behavior. Thus, if this is correct, one can definitely rule out a KT 
transition. However, as shown here, it is in fact not correct: in this Comment, we present 
simulation results for the fourth-order helicity modulus (T4) and from the numerical results 
we directly verify the discontinuous character of the helicity modulus at the upper phase 
transition of the six-state clock model. 

Reference |6j proposed a numerical method to identify the KT transition in the two- 
dimensional XY model based on a stability argument. That is, an external twist of mag- 
nitude A across the Xy-spin system gives an additional contribution to the free energy 
F so that F{A) > F{A = 0). For the system to be stable under small A, where one 
can expand the free energy density / (= F/N with A^ being the number of XY spins) as 
/(A) = (T) ^ + (T4) ^ + ■ ■ ■, the helicity modulus (T) must be non-negative. In the XY 
model, this quantity is zero above the critical temperature T^t and positive finite below T^t 
in the thermodynamic limit. The same stability consideration tells us that the fourth-order 
helicity modulus (T4) also must be nonnegative at any temperature T where (T) vanishes. 
Suppose that (T4) is finite and negative at the transition. Then (T) cannot approach zero 
continuously but must instead make a discontinuous jump to zero at the transition. Such 
a discontinuous jump is precisely the characteristic behavior of (T) at the critical tempera- 
ture for a KT transition. Unfortunately, it is notoriously difficult to verify the discontinuous 
character of the helicity modulus directly from numerical simulations because the precision 
of the simulations are restricted by the finite size of the system simulated. However, the 
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FIG. 1: (Color online) Temperature dependences of correlation functions, (a) the helicity modulus 
and (b) the fourth-order helicity modulus for the six-state clock model, (c) Amplitudes satisfying 
(T) = 2T/-K. The extrapolation yields (T) = 0.5742(3) at L — > oo. (d) Size dependence of the 
depth of (T4). Shown are L = 32, 64, 128, 256, and 512 from right to left both in (c) and (d). 

fourth-order modulus (T4) described above does not have this problem, as shown in Ref. [6|. 
We have here calculated the helicity modulus and the fourth-order modulus for the six- 
state clock model on L x L square lattices with the periodic boundary condition using the 
Wolff algorithm [7| [Figs. [U^a) and[T]^b)]. One clearly sees that (T4) approaches a negative 
finite value, which means that this system will exhibit a discontinuous jump in (T) at the 
transition temperature Tc according to the above argument. In addition, if it exhibits the 
universal jump J8|, the transition temperature and the jumping amplitude is related by 
(Y) = 2Tc/7i. Noting that the correlation length ^ scales as log^ ~ log-L ~ (T — Tc)^^^^ in 
the KT scenario, our extrapolation to L — t- 00 gives Tc = 0.9020(5) as shown in Fig. Wl^c). 
This agrees well with an independent estimation of T^ = 0.9008(6) in Ref. j4|, differing only 
by 0.3% at most. Recall also that Ref. [l| estimated T^ as 0.997(2) ignoring the KT scenario, 
which cannot be correct since one finds that the peak of (T4), which signals the criticality, 
is below this temperature for large enough L [Fig. [Il^b)]. These observations rule out the 
possibility that the upper transition can be of a non-KT continuous type as suggested in 
Ref. y. 



What could be the reason for the result in Ref. jl|? First of all, one notes that in Ref. [l| 
the same critical index v is found for both the upper and lower transitions, suggesting that 
the two transitions are identically continuous. However, this is established only for rather 
small lattice sizes, L < 28 used in the study, while we use lattice sizes up to L = 512 
here. Thus the results found in Ref. l| could be an artifact of the small lattice sizes. In 
particular, one may notice that for the generalized clock model studied in Ref. [9| the two 
separate transitions merge into one, when the cosine interaction is slightly distorted. At the 
merging point, the transitions can be described as one joint continuous transition. Therefore, 
one possibility could be that the two transitions for the usual six-state clock model are, for 
small enough lattice sizes, both strongly influenced by the critical behavior of the merging 
point. If this is the case, then the Fisher-zero approach should become consistent with a 
KT transition provided that the lattice sizes are large enough. 
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